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Abstract 

In this work, we report measurements of the height profile of transverse dunes 
in the coastal dune field known as "Lengois Maranhenses", northeastern Brazil. 
Our measurements show that transverse dunes with approximately the same height 
present a variable brink position relative to the crest, in contrast to the case of 
barchan dunes. Based on our field data, we present a relation for the dune spacing 
as a function of the crest-brink distances of transverse dunes. Furthermore, we 
compare the measurements with simulations of transverse dunes obtained from a 
two-dimensional dune model, where a phenomenological definition is introduced for 
the length of the separation streamlines on the lee side of closely spaced transverse 
dunes. We find that our model reproduces transverse dune fields with similar inter- 
dune distances and dune aspect ratios as measured in the field. 
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1 Introduction 



Dune morphology has been investigated over the last fifty years by scientists 
from different research areas like geologists, geographers, ph ysicists and even 
mathematicians, since the pioneering work of lBaenoldl (jl94ll) . Although a the- 
oretical explanation for the wide variety of observed dune forms is still lacking. 
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it is well known that the behaviour of the wind as well as sand availability are 
fundamental parameters for determining the final dune shape. On bedrocks 
and surfaces with not too much sand and under uni-directional winds, barchan 
dunes move separated from each other, without changing their shape and with 
a velocity inversely proportional to their height. These dunes have a windward 
side and a steep slip face where avalanches take place maintaining dune sta- 
bility and motion. As the amount of available sand increases, "barchanoids" 
are found, which look like oscillating chains of barchans joined at their horns, 
aligned perpendicularly to the wind direction. For higher sand supply, trans- 
verse dunes are observed. Transverse dunes cover 40% of all terrestrial sand 
seas, being common in the Northern Hemisphere, particularly in China, and 
dominating the dune fields on Mars. Normally transverse dunes have more ir- 
regular patterns and even hierarchies of smaller dunes can be f ound ori them . 
Transverse dune s have been studied in the field many times ( Wilsonl. 1972 ; 
Lancasteii . ll982L llQS.i iMuUiganl . Il988l : iBurkinshaw and RustL Il99.i IWigg, " 
but still most of these works are restricted to desert dunes. 



Dunes appear very commonly along coa sts fae sp et"aL . 1989*:*Nordstr om et al 



1990" Hesp . 2002,:. Kleinhans et al. . 2002^ Gongalves et al.l .i2003: Barbosa and Domingue^ . 
2004 ) . The sand that constitutes these dunes comes from the sea being there- 



after deposited on the beach. Once the grains are exposed to the air, they dry 
and some can be carried by the wind, initiating sand transport. Coastal dune 
development is a result of the complex interplay of several factors like the 
type of sediment, the sand influx, and the presence of vegetation. Important 
aspects for coastal dune formation are t he topograp hy of the area, the wave 
climate and, on a larger scale, sea level ( Hespl . 2002| ). With heights that can 
vary between a few meters to over 30 m, dunes represent not only tourist at- 
traction in coastal zones, but they also work as barriers against waves, act as 
sand reserve for beaches and protect the land behind them from erosion and 
salt intrusion. Understanding the morphology and dynamics of coastal dunes 
is thus of interest also for federal, state, and local governme nt units respons ible 
for the protection and management of coastal dune areas ( IBAMA , 2003| ). 



For strong unidirectional winds blowing from the sea onto the continents, 
transverse dunes appear along the coasts, moving landward. Barchanoids and 
transverse dunes extending over several kilometers form the characteristic pat- 
tern of t he biggest coastal dun e field in Brazil, known as "Lengois Maran- 



henses " (Goncalves et al. . 20031 ). Such scenario is typical of transgressive dune 



fields ( Hesp et al. . 1989| ). which appear on high energy beaches experiencing 
high littoral drift. In this work, we present field measurements performed 
on a sequence of transverse dunes situated at the beginning of the Lengois 
Maranhenses. The height profile of seven transverse dunes was measured and 
compared with numerical s imulations using a dune model in two dimensions 
fjSchwammle and Herrmann . .2004) , where in the present work a phenomeno- 
logical relation for the separation bubble has been introduced to the model. 
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based on the field data for interdune spacing as a function of the dune height 
and crest shape. 



2 Area of investigation 



The National Park of the Lengois Maranhenses, also known as the "Brazilian 
Sahara" , is located on the coastal area of the Maranhao State, in northeastern 
Brazil (Fig. 1), delimited by the coordinates S 02°19' and 02°45', and W 42°44' 
and 43°29'. It has a total area of 155 thousand hectares, and a coast of a length 
of 50 km, along which strong winds, reaching velocities of 70 km/h, transport 
sand from the Atlantic Ocean onto the continent. The wind rose in Fig. 2 
shows that winds blow mainly from the East. 

The area of the "Lengois" (or 'sheets' in Portuguese) is mainly characterized 
by the presence of barchanoids and transverse sand dunes separated in the 
rainy season by lakes and lagoons (Fig. 3). The area has a semi-humid tropical 
climate, with little vegetation. The local temperatures are typically 30 to 38°C, 
with annual rain between 1500 mm and 2000 mm, which is an amount quite 
different from the mean rainfall in desert areas (in general less than 250 mm). 
The rainfall distribution, however, is very concentrated, as we can see in Fig. 
4. In the first months of the year, high rainfall makes sand transport difficult, 
thus implying low dune mobility. A large fraction of the interdune lakes formed 
in this period disappear in the dry season. The landscape of the dune field 
appears to change continuously, indeed dune mobility implies that the lakes 
often reappear in different places with different contours. The lowest rainfall 
indices in the studied area refer to the period between August and November. 
The field work has been carried out from 23 to 29 of September of 2003, in 
the middle of the dry season. 

The transverse dune field studied is shown in Fig. 5. It is situated in the 
vicinity of Atins, a small village in the Maranhao State (MA), just 1 km from 
the sea. These dunes have a height of typically 10 m and are closely spaced, 
separated by a few meters only, on a dense sand sheet in the beginning of 
the Lengois Maranhenses. They may be classified as "transgressive dunes", 
which are defined in most of the cases as dunes that develop on sand sheets 
f)HesD et al.l . ll989l) . 



Near to the transverse field, several small barchans and barchanoids are found, 
with heights ~ 50 cm — 1 m. These smaller dunes emerge on a fiat surface 
covered with less sand, extending from the beach up to the transverse dune 
field, as shown in Fig. 5(b). The size difference between the two neighbouring 
dune groups is clearly visible. The dune field studied here represents the gen- 
esis of the Lengois, where the characteristic transverse profile originates. The 
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Fig. 1. Map of the LeriQois Maranhenses region. The studied transverse dune field 
is located near the village of Atins, MA. 
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Fig. 2. Wind rose for the region of the Lengois Maranhenses corresponding to the 
period from January to December of 2003. We see that the wind blows mainly from 
the East. The velocity of the wind ("u") is shown in units of m/s. 




Fig. 3. Typical landscape of the Lengois Maranhenses. Chains of barchanoids are 
separated by interdune lakes that emerge in the wet seasons. Photo: Morals Brito. 

dark areas in Fig. 5(b) correspond to low vegetation, called "mangues". The 
coastal zone near Atins is also characterized by a high concentration of plants 
developed on parabolic dunes of different sizes. Vegetation appears due to hu- 
midity from the ocean and interdune lakes, and grows during the wet seasons, 
appearing also, although more rarely, on some parts of transverse dunes. 



In the next section we present the results of our measurements. In Section 4, we 
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Fig. 4. Rainfall throughout the year. Vertical axis is in mm month 



descri be the dune model originally introduced in ISchwammle and Herrmann 
( 2004^ and propose a modification of the separation bubble to simulate closely 
spaced transverse dunes, based on the data of our measurements of dune 
spacing and crest-brink distances of the dunes. Comparison of the results of 
our simulations with the measured profile and discussions are presented in 
Section 5. Finally, in Section 6 we present our conclusions. 



3 Measurements 



A tachimeter was used to measure the height variations at 77 different points 
along the surface of the dunes, which define a profile with a length of almost 
720 m, as shown in Fig. 6. To interpolate the entire profile of the transverse 
dune field, we used the relative distances and height differences between the 
measured points and a fixed angle 6c = 34° for the slip face for each dune. 
The reference height /i = is set at the foot of the slip face of dune 3, whereas 
X = corresponds to the brink of dune 1, for which we did not measure the 
whole profile. Our measurements have error estimates of 0.1 m for the height 
coordinate and 1 m for the horizontal axis. The measured transverse dunes 
have heights of typically 7 — 10m and a crest-to-crest distance, Dec, of 90 — 130 
m, as shown in Table 1. The sixth column of this table shows the aspect ratio r 
of these dunes, which is obtained by dividing their heights at the crest, i^crcst; 
by the corresponding values of the length of the windward side, Lq, and we 
find that r = 0.08 ± 0.02 is consistent for all dunes. The quantities listed in 
Table 1 are defined in Fig. 7. 

It is surprising that although the dunes have more or less the same height, the 
position of the brink with respect to the crest, dc, varies strongly. This is in 
strong contrast with the situation of single barchans, for which it has been ob- 
served from field measurements in sout hern Morocco that the br ink gets closer 
to the crest for higher barchan dunes ( Sauermann et al. . 2000l ). In addition. 
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Fig. 5. In (a) we show a Landsat image (1984) of the Lengois Maranhenses, where 
the studied transverse dune field is indicated by the arrow. In (b) we show an image 
of the studied field, from Embrapa Monitoramento por Satelite, 1999. 

differences between small and large barchans h ave been recently report ed by 
Hesp and Hasting^ ( 1998| ). On the other hand, Long and SharpI ( 1964 ) have 
more generally identified two kinds of barchan dunes: the "slim" and "fat" 
barchans, where the latter dunes may have been the more "morphologically 
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Fig. 6. Height profile of the transverse dune field studied. 



Dune 


-ff crest (m) 


-ff brink (m) 


dc(m) 


Lo(m) 


r — H crest /Lq 


Dcc{i,i + l)(m) 


1 




6.6 ±0.1 










2 


9.6 ±0.1 


6.6 ±0.1 


32 ± 1 


123 ± 1 


0.078 ± 0.005 


132 ± 1 


3 


8.8 ±0.1 


8.6 ±0.1 


7± 1 


96 ± 1 


0.092 ± 0.005 


94 ± 1 


4 


7.0 ±0.1 


7.0 ±0.1 


0± 1 


65 ± 1 


0.108 ±0.005 


95 ± 1 


5 


6.9 ±0.1 


4.9 ±0.1 


24 ± 1 


95 ± 1 


0.073 ± 0.005 


119 ± 1 


6 


9.7 ±0.1 


9.7 ±0.1 


0± 1 


88 ± 1 


0.110 ±0.005 


124 ± 1 


7 


7.9 ±0.1 


4.3 ±0.1 


35 ± 1 


127 ± 1 


0.062 ± 0.005 





Table 1 

For each dune in the profile of Fig. 6 we identify in columns 2, 3, 4, 5 and 6, 
respectively: The dune height at the crest, .ffcrest, the height at the brink, -ffbrinkj 
the crest-brink distance dc, the windward side length Lq and the aspect ratio r. The 
seventh column shows the crest-to-crest distance Dcc{i, i + between the i—th and 
(i ± 1)— th dunes, -ffcresti -f^brinki and Lq are shown in Fig. 7 for dune 5. 




Fig. 7. Profile of dune 5 in Fig. 6 with definition of variables Merest, -f^brinki dc 
Lq in Table 1. Heights {H) are measured relative to the foot ot the slip face. 



and 



complex and areally larger dune masses" form ed by coalescen ce and fusion, 
and behave in a different way as the slim ones. iBagnoldl (|l941l ) regarded this 
phenomenon as the product of the lag in-between changes of wind conditions 
and the resulting changes in sand movement. 



The upper inset of Fig. 8 shows the inclination 9^ of the windward side (circles) 
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Fig. 8. Difference between dune heights at the crest and at the brink, as a function 
of the distance dc between the referred two points for the measured dunes. For 
dunes 4 and 6 the crest coincides with the brink, while the corresponding data for 
dune 1 were not recorded. The observed quadratic relation is a consequence of the 
parabolic shape of the dune at the crest. The upper inset shows the aspect ratio 
r = i^crest/-^o = 0-08 it 0.02 (triangles) and the average windward side inclination 
tan 6^ = Hcrest/{Lo — dc) = 0.10 ± 0.02 (circles), 0^ ~ 6°. The lower inset shows 
the crest-brink distance versus the height of the crest. 

together with the aspect ratio r (triangles) of the measured transverse dunes. 
Here 9^ is defined through the relation tan^^ = i7crest/(^o ~ ^^c); where dc 
is the horizontal distance measured from the crest to the brink of the dune. 
From the value tan^^„ = 0.10 ± 0.02, we get 6^ ^ 6°. The main plot of Fig. 8 
shows that the difference between the heights at the crest and at the brink is 
essentially quadratic in dc- Moreover, Fig. 9 shows that the transverse dunes 
present a nearly parabohc shape at their crest. In this figure, we lump the 
profiles of dunes 2 — 6 together, and we can see that the windward sides of all 
dunes at their crests superimpose in a single parabolic curve with maximum 
at point C. The vertical (horizontal) axis of the main plot gives the profile 
H (horizontal distance x) measured relative to point C for all dunes. The 
empty circles in the inset of Fig. 9 represent the sampled profiles for all dunes 
close to point C, as indicated by the dotted lines, while the continuous line 
corresponds to the curve y — —kix"^, with ki = 0.0025. 



For dune 3, we performed correlated measurements of the wind velocity v and 
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Fig. 9. The main plot shows the height profile of dunes 2 — 6, where the dunes are 
identified by numbers at their slip faces, and are represented by different symbols. 
The windward sides of these dunes are shown superposed to each other in such a 
way to form a parabola with maximum at C(0, 0). This behaviour is more clearly 
shown in the inset, where the circles represent the height profile of all dunes of the 
main plot, and are fitted by a parabola (continuous line). The symbol "~" over H 
and X indicates that these variables are measured, both in the inset and main plot, 
relative to point C for each dune. The dotted lines in the inset and main plot are 
guides to the eye meant to localize the maximum C of the parabola. 

the sand flux q on two different points along its windward side. This dune 
has a height -ffcrest = 8.8 m measured from the foot of its slip face up to 
the crest, and the length of its windward side is approximately 100 m. We 
placed an anemometer at a height of Im over the ground, and the registered 
wind velocity was averaged over time intervals of 10 min. An average wind 
velocity of (i;a) = 7.5 ± 1.2m/s was registered at the brink of the dune (point 
"A"). Then, we moved the anemometer to a point at about the middle of the 
windward side (point "B"), and the average wind velocity registered there was 
{vb) = 5.4±0.5m/s. A second anemometer, placed at one point outside of the 
dune fleld, at about 10 m from the transverse dunes on a sand sheet where 
a fleld of small barchans developed close to the ocean, registered an average 
wind velocity of (wr) = 5.3 ± 0.7 m/s, also at a height of 1 m. We did not 
measure the wind velocity at different heights over the ground, and thus we 
could not estimate the shear velocity u^q and the aerodynamic roughness zq 
of the surface. This quantities can be obtained from the logarithmic proflle of 
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Fig. 10. Measuring wind velocity and sand flux using an anemometer and the three 
traps at the different points of dune 3 in Fig. 6. 

the wind velocity v{z) (jPve and Tsoail . Il99(t : 



v{z) 



— m — . 

K Zq 



where z is the height over the ground and k = 0.4 is the von Karman constant, 
and M* is used to define the shear stress of the wind: 



r 



(2) 



where pair is the density of the air. If we take, for example, a roughness length 
zq = 2.5 X 10~^ m, we then obtain the undisturbed shear velocity over the plane 
~ 0.36 m/s, using v{z = Im) = 5.35 m/s. It is interesting to note that 
these were essentially the same values of m*o and zq found i n a recent work on 



the dune field of Jericoacoara, also in northeastern Brazil ( Sauermann et al 



2003| ). Since we could not measure accurately the values of these quantities, 
we use in our simulations (see Section 5) estimated values of u^o close to 0.4 
m/s and zq between 0.8 — 3 mm. 

At the same positions of the dune where the anemometer was placed, as shown 
in Fig. 10, we fixed 3 cylindrical traps to measure the corresponding sand fiux. 
Each trap consists of a cylinder of 1 m height and 4.8 cm diameter, and has an 
open vertical slot of 1 cm width. The traps were turned with their slit facing 
the wind and their bottom buried in the sand (see Fig. 10). The back of the 
traps has an opening of 2 cm width covered by a permeable fabric with pores 
smaller than the grain diameter. Initially, the traps were placed with their 
back showing in the wind direction. Then, they were simultaneously turned 
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to align their op enings with the wind dir ection and sand was collected, in a 
similar way as in Sauermann et al. ( 2003| ). From the mass m of the collected 
sand, the width w of the opening and the time T necessary to fill the traps, we 
calculated the sand flux q = m/(Tw), for which we use m = irr'^hc x Pquartz, 
with r = 2.4 cm and Pquartz = 2650 kg/m^. For each point two traps were used 
which captured grains travelling up to a height he ^ 17 cm and another one at 
a height of 1.5 cm above the ground. The mean values (q) obtained from two 
runs using these traps are, for points "A" (brink) and "B" (middle of windward 
side) respectively, (^a) = 0.023±0.005 kg/m-s and (gs) = 0.009±0.001 kg/m-s, 
where the error estimates are defined by the maximum and mininum values 
registered for q at each point, as in the case of the measurements of the wind 
velocity. A rigorous treatment of our measurements of sand flux would require 
reference to a recent work on sampling efficiency of a vertical array of aeolian 
sand traps ( Li and Ni . 2003[ l. Since we have only two measurement points for 
wind and flux, we report here the measured values but in this work we do not 
show comparison with model predictions for these data. We present typical 
simulations of the transverse profile and compare our results with the observed 
dune shape. 



4 A two-dimensional model for transverse dunes 



A two-dimensional mod el including the m ain processes of dune morphology 
was recently intr oduced (iKrov et a l.l '2002|) a nd ex tended to study transverse 
dunes in a field ( Schwammle and Herrmannl . 120041 ). From this model, the air 
shear stress r and the saltation fiux q at the ground are numerically calculated 
in each iteration, as well as the fiux due to avalanches at the slip face of 
the dunes. The shear stress perturbat ion is calculated in the two-dimensional 
Fourier space using the algorithm of IWeng et al.l (|l991[ ) for the components 



Tr and r„ 



\k\U^{l) 



1 + 



2 ln(L|fc^|) +47 + 1 + isign(fc^)7r' 
In (l/zo) 



(3) 



and 



. , s _ 2h{kx,ky)kxky 



where the axis x (y) points parallel (perpendicular) to the wind direction, k^ 
and ky are wave numbers, \k\ = ^k^ + ky, 7 = 0.577216 (Euler's constant), zq 
is the roughness length and L is the characteristic length for a hill ({Hunt et al.[ 
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1988| ) , and is defined as the length measured from the half height of the wind- 



ward side to the crest. U{1) = v{l)/v{hra) is the undisturbed logarithmic profile 
(1) calculated at height /, which is given by 



and normalized by the velocity at a reference height = Lj y'logL/zo, 
which sepa. r ates t he middle and upper flow layers defined in the model by 
Hunt et al. ( 1988| ). It has been verified from simulations that an open sys- 
tem of transverse dunes seems to reach lateral invariance under unidirec- 

Although this observation is 



200J 



tional wind (jSchwammle and Herrmann , 
not strictly valid in the field due to the multi-directional nature of secondary 
airflow patterns ( Walker and Nickling . I20n2[ ). the assumption yields to im- 
portant advantages in transverse dune field calculations. Thus, for modelling 
transverse dune fields, as opp osed to the situation of barchan dune modelling 
JSchwammIe and Herrman;i . B , we neglect the ^/-coordinate and consider 
height profile variations in the x direction only, since much less computational 
time is consumed in this case, and larger dune fields can be simulated. In this 
way, we use the shear stress perturbation evaluat ed for a slice of a transverse 



dune translationally invariant in the y direction ( Krov et al. . 2002f l: 



r X 



\tx J X — E 



(6) 



where h' = dh/dx is the height derivative in x direction, and the constants 
fi{zo/L) and f2{zo/L) depend logarithmically on the ratio zq/L. The shear 
stress is then given by: 

r = r(x) = ro[l + f(x)], (7) 



where Tq is the undisturbed air shear stress over a flat ground. Of course, 
the fluctuations in the y— direction that can be observed in Fig. 5(b) may be 
of further importance for the dune dynamics, but are not considered in this 
approach. The sand flux q is calculated as: 

dp ^ d 




where p and u are, respectively, the mean density and the mean velocity of the 
saltating grains, and are both integrated in vertical direction. The saturated 
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density at the steady state, ps, and the characteristic time that defines the 
transients of the fiux are given by: 



Ps 



2a 

-ir-n) 



and T, 



2au 



9 l{r-Tt) 



(9) 



where rt is the threshold shear stress for saltation transport (jBaenoldl 119411 ). 
and a is a model parameter that represents an effective restitution coefficient 
for the grain-bed interaction. The vertical component of the initial velocity of 
the ejected grains after the impact of saltating grains onto the bed ("splash") 
is equal to the average gain in horizontal velocit y of the impacting grain, after 
its acceleration by the wind, multiplied by a ()Sauermann et al.l . 



the other han d, 7 is the second model parameter, to which the erosion rate is 



proportional ()Sauermann et al.L 120011 1. 



The steady state is assumed to be reached instantaneously, since it corresponds 
to a time scale several orders of magnitude smaller than the time scale of the 
surface evolution. Thus, time-dependent terms are neglected. Furthermore, the 
velocity u of the grains in the saltation layer is modelled through the following 
equation: 



3 Pair 1 . 

jCd -{VeS 

^ Pquartz 



u)\VcS - U\ 



9_ 

2a 



0, 



(10) 



where: 




:ii) 



and u^t = ^Tt/pair- In Eq. (11) Vcs is an effective velocity calculated at a ref- 



sand 



erence height Zi within the saltation layer between the roughness length z\ 
associated with the grain diameter and the mean sa ltation height z-m- The value 
of 

^sand jg smaller than the aerodynamic roughness (jPve and Tsoan. Il990> ). and 
is taken as of the order of d/10, where d is the grain diameter. The velocity u is 
calculated from t>efr using eq. (10). Typical values of the model parameters zi, 
a and 7 defined above have been determined from comparisons of the model 
with microscopic simulatio ns and wind tunne l meas urements of sand fiux and 
saturation time transients (|Sauermann et al.Ll2001[ ). We refer to these results 
and use the same value s for the constants and microscopic model parame- 
ters a s in previous work ( Sauermann et al. . 200lt Schwammle and Herrmann . 
20041 ): g = 9.81 m/s^. 



K 



0.4, 



Pa 



= 1.225 kg/m3, pq,,rt, = 2650 kg/m^, 
Zl = 0.005 m, Zm = 0.04 m, a = 0.35, 7 = 0.4, ~ 0.3 m/s, d = 250 pm and 
Cd = 3. It has to be emphasized that these values have not been obtained from 
measurements on the field studied in the present work. Of course, a detailed 



14 



quantitative description of the sand flux and wind velocity over the transverse 
dunes investigated here would require more field data. Furthermore, the shape 
and dynamics of dunes may be strongly influenced by the sorting and size dis- 
tribution of grains, which is in general heterogeneous rat her than well defined 
at one single size, and is also related to the dune's age (Besler, 2002| ). In our 
model, such a feature is not taken into account. A recent study has shown that 
the mean grain diameter of the san d in the area in vestigated in the present 
work lies between 177 and 354 /im fliBAMAi bnnat ). In simulations, a grain 
diameter of d = 250 /im is typically used, around which the grain size d i stribu - 
tion of sand dunes appears sharply peaked as shown by Pve and Tsoad (Il990f) 



thus definin g a threshold shear velocity of about 0.3 m/s ()Sauermann et al 



2OOILI2OO3I ). 



For the avalanches at the slip face, the BCRE model (|Bouchaud et al.L 11994 ) 
is used to calculate the new surface after the relaxation by sand transport 
downwind: 



at \ ox 



tan O 



:i2) 



OR d ^ ^ ^(dh 
— + T^Rua = CaR{ ^ - tane 
at ox \ox . 



(13) 



where h is the height of the non-moving sand bed, R is the height of the 
moving layer, Ca is a model parameter and Ua is the velocity of the grains in 
the moving layer. Avalanches are considered to occur instantaneously, since 
the time scale associated wit h them is negligible if compared w ith the time 
scale of the surface evolution (|Schwammle and HerrmannL 12004^ . 



The evolution of the dune surface is determined by erosion and deposition of 
sand, and its profile can be obtained from the mass conservation equation: 



dtp + 



9$ 

dx 



(14) 



where the sand density, p, and the sand flux per time unit and area, $, are 
integrated over the vertical coordinate assuming that the dune has a constant 
density of Psand: 



1 

Psand 



pdz, q= ^dz, (15) 
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where psand = 0.62pquartz, leading to the following equation for the surface 
height h: 

dt Psand dx ' 



Equation (16) is used to determine the characteristic time scale of the model, 
typically 3 — 5 hours for every iteration. Simulation steps may be summarized 
as follows: (i) the shear stress over the initial surface (a sand bed with small 
fluctuations as described below) is calculated using Eqs. (6) and (7); (ii) from 
the shear stress, the sand flux is calculated with Eq. (8); (iii) the change in 
the surface height is computed from mass conservation (Eq. (16)) using the 
calculated sand flux; and (iv) if the inclination of the surface gets larger than 
34°, avalanches occur and a slip face is developed. Steps (i) — (iv) are iter- 
atively computed, and in case of transver se dunes, the evolution of the dun e 
field depends on the boundary conditions ( Schwarnmle and Herrmannl. 12004 ) . 



For pe riodic boundary conditions, it was shown in lSchwammle and Herrmann 



f)2004^ that the dunes increase continuously in height though they keep their 
average aspect ratio, thus comparison with real profiles are made taking sim- 
ulation snapshots of similar dune size as observed in the field. 

For a dune with slip face, flow separation occurs at the brink, which repre- 
sents a discontinuity of the surface. The flow is divided into two parts by a 
streamline connecting the brink with the ground at the reattachment point. 
These streamlines define the separation bubble, inside which eddies occur and 
the flow is often re-circulating, and may present a complex multi-directional 
pattern ( Walker and Nicklind . 20021 ) , whereas the averaged flow outside is lam- 



inar. In the area of the separation bubble, we set the shear stress equal to zero, 
and in our model no sand transport occurs in the bubble. The length of the 
separation bubble is thus of major importance for interdune spacing. 

We observed in the field that the net interdune sand fiux in the wind direction 
was zero from the foot of the upwind slip face up to the beginning of the 
windward side of the downwind dune. If we measure the horizontal distance 
between the brink of dune i and the point at which the windward side of dune 
i + 1 begins, i.e. Xj+i — Xbrink,?. for the dunes in the measured field, we find that 
this quantity is between 1.5 and 4 times the height of dune i at its brink, -ffbrink.i 
(see inset of Fig. 11). These interdune distances can not be reproduced with 



the or iginal transverse dune model proposed in ISchwammle and Herrmann 



(|2004[) . where the separation bubble was defined in such a way as to have 



lengths between 6 and 10 times the dune height. Here we have to remark 
that the separation bubble defined in our model is the region with zero net 
fiux at the lee side of the upwind dune. As mentioned above, the actual wind 
behaviour and, thus, the fiux at the lee si de have been found to presen t complex 
patterns before the reattachment point ( Walker and Nickline . lioO^ . a feature 
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Fig. 11. Main plot: Distance between the slip face foot of dune i and the windward 
side foot of dune z + 1 as a function of the distance measured between the brink 
and the maximum of the parabola according to Fig. 9. As we can see, the length 
of the separation bubble appears to decrease linearly with Abrink-crest,i- The inset 
shows the distance Xj+i — Xbrink,i between the foot of the windward side of dune i + 1 
and the brink position of dune i as a function of the height of dune i at the brink, 
-ffbrink,i- The dotted line in the inset represents the curve y = 4 x (iJbrink.i — 2.5), 
which is a guide to the eye meant to show that interdune distances are below four 
times the dune height at the brink. 

that is not captured in our model. On the other hand, as we can see in Fig. 11, 
interdune spacing also depends strongly on the position of the brink relative 
to the crest. Here, we plot the horizontal distances between the dunes, A^ j+i, 
measured from the foot of the slip face of dune i to the beginning of the 
windward side of dune i + 1, as a function of the distance between the brink 
of dune i and the maximum of the parabola in Fig. 9. In other words, dunes 
that have the position of the brink at the left (right) of the maximum of 
the parabola (point C) in Fig. 9 are associated with a proportionally larger 
(shorter) interdune spacing relative to the corresponding downwind dune. 



IParsons et al. ( 2004 ) have found lengths of 3 to Vbh for the separation zone 
using CFD simulations, where h is the height of the transverse dune. However, 
in this work, they only studied sharp crested transverse dunes, where the 
crest coincided with the brink. Simulations ai id measurements of surface shear 
stress over transverse dunes in a wind tunnel ( Walker and Nickline . 2003| ) have 
reported similar reattachment lengths for such sharp crested dunes. In a recent 
work, the profile of the wind over the transverse dunes studied here has been 
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calculated using FLUENT ( Herrmann et al. , 200^ . It was shown in this work 
that there is in fact a strong dependence of the length of the separation bubble 
with the shape of the dune at its crest, where dunes with a large crest-brink 
distance, i.e. very r ound dunes, almost do n o t show recirculating flow at its 



lee side. Moreover, iMcKenna Neuman et al. have noticed differences 



in the behaviour of the streamlines, flow expansion and deposition for sharp 
crested dunes in relation to the rounded ones. 



The observations above lead to the conclusion that simply considering the 
length of the streamlines as proportional to the height at the brink is not 
reasonable for modelling transverse dunes, and that the observed dependence 
of the interdune spacing on the crest-brink distance of the dunes has to be 
taken into account. Since in our model the net sand flux is set to zero in the 
separation bubble, the observed interdune distances define an upper bound 
for the reattachment length. If we choose a larger length for the bubble in 
our model, simulation of closely spaced transverse dunes as the ones mea- 
sured here is not possible (jSchwammle and Herrmannl . |200^ . The linear rela- 
tion between the interdune spacing and the distance Abrink-crest between the 
crest and the brink, which is shown in Fig. 11, suggests that, if the horizon- 
tal distance between brink and crest is described by a parabola of the form 
ifcrest — -f^brink = ^1 A^rink-crest ' then in a first approximation we could define 
the distance from the crest to the reattachment point to follow also a quadratic 
relation -ffcrest = fc2(a^reat — a; crest This relation is then used to determine the 



point of reattachment x^eat of the separation bubble, Xreat = y H^est 1^2+ x crest- 
In order to give an insight into the meaning of the parameter we no- 
tice that for a dune with height 7 m at its crest, a value of k2 = 0.02 gives 

a^reat " a;crest = 18 m, whilc k2 = 0.01 leads to Xreat - a;crest = 26 m. 



In our previous work, the streamlines of the separation bubble were fitted by a 
third order polynomial, /(x) = a^x^ + 02^^ -|- aix + Oq, whose parameters were 
calculated from the continuity conditions of the profile h{x) at the brink of the 
dune and at the reattachment point, as well as of dh/ dx at the brink. Moreover, 
the third order polynomial was calculated i n such a manner that the separatio n 
bubble doesn't descend steeper than 14° (jSchwammle and Herrmannl . l2004 ) . 
This leads to large values of the horizontal distance from the brink up to the 
reattachment point, Xreat — a;brink being typically between 6 and 10 times the 
dune height. Using higher order polynomials to take into account the relation 
observed in Fig. 11, on the other hand, makes the streamlines descend too 
fast close to the brink at the lee side, making the curvature of the separation 
bubble different from that of the dune surface at the brink, thus giving rise to 
incon sistencies in the ca lculation of the shear stress with Eq. (6) for smooth 



hills (jWeng et al.l . ll99lh . and consequently of the sand flux. 
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Fig. 12. Profile h{x) of measured dune 3 with the modified separation bubble s{x). 
The 3rd order polynomial separation bubble f{x) is shown for comparison. 

Therefore we define a modified separation bubble as = /(x) for Xbrink < 
X < Xcnt and 



S{x) = [/(x) - /l(Xreat)] X {1 " tanh [{X - R)/\c]}/2 + /l(Xreat), 



(17) 



for Xcut < X < Xreat, where X^nt = a^brink + C{x reat ^^brink); R — (^^cut ~l" 

a;reat)/2, Ac determines the slope of the streamlines at point R, and C is a 
parameter between and 1, which determines the point at which the original 
polynomial is multiplied by the function tanh. The lengths involved in the 
modified separation bubble are schematically shown in Fig. 12. 



Eq. (17) gives a shorter separation bubble from the original function f{x). 
It keeps the slope of the original separation bubble close to the brink, and 
presents a steeper descent to the reattachment point. Thus, E g. (17) is a simple 
appro ximation of the original separation bubble introduced in lSauermann et al 
( 20011 1 to study closely spaced transverse dunes, based on the phenomenolog- 
ical observations reported in the present work. We set ^2 = 0.014, which gives 
a distance of 25 m from the crest of dune 3 up to the reattachment point at 
its lee side, and Ac and ( are set to 0.1 and 0.9 respectively. 



Studies of the separation bubble of isolated transverse dunes have shown that 
the separation streamlines may be well described by an elliptical curve, but still 
separation distances only larger than 4 times the dune height were found, and 
no mention has been made to closely spaced transverse dunes in a field. Besides 
this, we report in the p resent work larger crest -brink distances than the ones 
of the dunes studied in Schatz and HerrmannI (j^05). Simulations using our 
model with the separation bubble lengths calculated in ISchatz and Herrmann 
(200i) would yield much more widely spaced transverse dunes than the ones 
investigated here. An extension of this work for transverse dune fields with 
dunes of different heights and crest-brink distances could yield more precise 
knowledge of the separation streamlines of closely spaced dunes. 
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In the next section, we present simulations of transverse dune fields in two 
dimensions and compare the results with the measured profile in Fig. 6. 



5 Simulations of transverse dune fields in two dimensions 



The numerical model described above has been used in Schwammle and Herrmann 



f)2004^ to study transverse dune formation on a sand bed completely filled with 
sand. As mentioned before, such dune field s developed on san d sheets have 



19891 ). For periodic 



been classified as transgressive dune fields (jHesp et al. 
boundary conditions, we choose at t = a sand bed of depth 10 m randomly 
disturbed by small Gaussian hills. We use the separation bubble introduced 
in the last section, which is able to give shorter distances between dunes than 
in the case ISchwammle and HerrmannI (j200J) . 



Figs. 13(a) to (f ) show simulation snapshots of a typical transverse dune field of 
2 km length obtained with our two-dimensional model. The model parameters 
are those defined in Section 5. The values of the shear velocity u^q and its 
threshold value u^,t are respectively 0.44 and 0.32 m/s. For comparison, we 
show in Fig. 13(g) the measured profile. The mean slope of the windward side 
in Fig. 13(f) is between 0.065 and 0.09, which are values close to the windward 
side inclination found in the measured field (0.08 — 0.12). 



Furthermore, we see in Fig. 13 that the distance between crest and brink for the 
simulated transverse dunes of similar heights varies strongly as in the case of 
the measured dunes. This striking observation appears to be a consequence of 
the interaction between transverse dunes in a field. Barchans behave more like 
isolated dunes and ex change less sand than a field of clos ely spaced transverse 
dunes. Simulations in lSchwammle and HerrmannI (|2004^ have shown that the 
brink position of transverse dunes is in general closer to the crest if the dune is 
downwind of another one. Such sharp crested dunes with "triangular" cross- 
section are found in our simulations during the evolution of the field, as shown 
in Fig. 13(c). However, they appear only when interacting with another dune 
like solitary waves, a feature that can be ilustrated comparing Figs. 13(c) and 
(d). In order to investigate dune interaction in more detail, transverse dunes 
have to be still systematically st udied as done recently for barchan dunes 
ISchwammle and Herrmann . 2003[ ). 



1 



The typical dune shape found in our simulations is the " rounde d " one , i.e. 
where the brink is clearly distinguishable from the crest. iTsoai (|198,^ has 
shown that the 'rounded' dunes are 'more stable'' than the 'triangular' ones, 
since the sharp-crested shape experiences strong flux at the brink. We have 
found such peaks of sand flux when trying to simulate the evolution of trian- 
gular transverse dunes with our model, and we have seen that they lose this 
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Fig. 13. From (a) to (f) we show snapshots of the simulation of a transverse dune 
field obtained from our two-dimensional model. Model parameters are mentioned 
in Section 4. The shear velocity u*o and its threshold arc respectively 0.44m/s 
and 0.32m/s. For comparison, the measured profile is shown in (g). 
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shape after some iterations, becoming 'round'. On the other hand, we could 
not find dunes with the 'triangular' shape in our simulations of dune fields, 
but those interacting like solitary waves, as mentioned above. However, we do 
not conclude from this that the measured sharp crested dunes found in Fig. 
6 should be interacting like solitary waves with their neighbours. We remark 
that our model assumes lateral invariance, from which the profile variations 
found for these dunes in the y direction, as we see in Fig. 5, are not included 
and, consequently, variations in the lateral shear stress are not considered, 
which could infiuence the dune shape at the brink. In spite of this limitation, 
the variable crest-brink distances found in the field have been captured by 
the model. Also, the separation bubble in our model is the region of zero fiux 
downwind, which implies that the observed interdune distances provide an up- 
per bound of the separation length. To use a longer bubble, our model could 
be impr oved in order to capture the complexity of the wind and fio^y in th e 
bubble (|Frank and Kocurekl . Il99fil : IWiggsl . bOOlt IWalker and Nickhnd . 120021 ) . 
Possibly, in this case, sharp crested dunes could be found in simulations of 
closely spaced transverse dunes. 



The shape of coastal dunes like the ones investigated here may be infiuenced 
by standing water and vegetation in the interdune area, which could be rele- 
vant to the lee face steepness and crest-brink separation distances. Recently, 
an extension of our model to simulate barchan dunes with vegetation has suc- 
cessfuU y reproduced the ob servation that th ese dunes t ransform into parabolic 



dunes ( Duran et al. . 2005[ ). as suggested in Tsoar and Blumberg (2002,). The 



study of coastal transverse dune field formation in the presence of vegetation 
and lagoons is the subject of further work. 



6 Conclusions 



We have performed field measurements of the height profile of transverse dunes 
in the Lengois Maranhenses, northeastern Brazil, and compared the results 
with the predictions of a two-dimensional dune model recently introduced. 
The model separation bubble at the lee side of the dune was adapted in order 
to match the observation that the interdune spacing in the measured field is 
less than 4 times the dune height at the brink. 

The shape of transverse dunes was found to vary with their height in a differ- 
ent way than in the case of barchan dunes. While larger barchans present their 
crest closer to their brink than smaller ones, the crest-brink distance was found 
to be independent of the height for transverse dunes. Moreover, we found that 
the windward side of all dunes could be superposed in such a way to form a 
nearly parabolic shape with common maximum point (see Fig. 9), indepen- 
dently of the relative position of the brink for each particular dune. Based on 
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the relation observed for the interdune distances as a function of the brink 
position, we defined the reattachment length at the dune lee side to simulate 
transverse dune fields similar to the measured one. Our model successfully re- 
produces the aspect ratio of the measured dunes, interdune distances and the 
appearence of variable crest-brink distance, assuming lateral invariance of the 
profile, which allows considerable simplifications of the model and simulations 
in two dimensions. In order to consider the oscillations in the profile in the 
y direction, as Fig. 5(b) shows, three dime nsional simulations are req uired to 
take into account the lateral shear stress fjSchwammle and Herrmann . 2005) . 
which could of course influence the dune shape for instance at the brink. The 
simulated dunes compared with the measured ones appeared 'rounded', and 
we have discussed possible explanations for this. 

Finally, a small field in the beginning of a growing structure of transverse dunes 
was measured which characterizes the pattern of the Lengois Maranhenses, 
as illustrated in Fig. 5(a). We have simulated a transve rse dune fie l d tha t 
evolved on a sand sheet, thus a "transgressive" dune field ( Hesp et al. . 19891 ). 



One aspect our model did not take into account is the proximity of the water 
line to the surface due to the rainfall in this region and the eventual role of 
standing water or vegetation, which may exert a major influence on the dune 
shape, although in the dry seasons water is practically absent from the sand 
sheet where the measured transverse dunes were found. Comparison of our 
field data with further field measurements on transverse dunes in the Lengois 
Maranhenses as well as in other areas with different climate and geological 
conditions could, thus, be of importance for better understanding formation 
and evolution of transverse dune fields. 
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